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Abstract. The isometric embedding problem is a fundamental 
problem in differential geometry. A longstanding problem is con- 
sidered in this paper to characterize intrinsic metrics on a two- 
dimensional Riemannian manifold which can be realized as iso- 
metric immersions into the three-dimensional Euclidean space. A 
remarkable connection between gas dynamics and differential ge- 
ometry is discussed. It is shown how the fluid dynamics can be used 
to formulate a geometry problem. The equations of gas dynamics 
are first reviewed. Then the formulation using the fluid dynamic 
variables in conservation laws of gas dynamics is presented for the 
isometric embedding problem in differential geometry. 



1. Introduction 

We are concerned with isometric embeddings or immersions (i.e., re- 
alizations) of two-dimensional Riemannian manifolds in the Euclidean 
space M^. A classical question in differential geometry is whether one 
can isometrically embed a Riemannian manifold {Ai"',g) into for 

large enough. Nash [12] indicated that any smooth compact man- 
ifold (A^", g) can always be isometrically embedded into M3s„+4ra 
Sn = n{n — l)/2. This important paper lays the foundation for the 
development of geometric analysis in the second half of the 20th cen- 
tury. Gromov [5j| proved that one can embed any (Ai", g) even into 
]^s„-i-2n+3^ Then a further natural question is to find the smallest di- 
mension N{n) for the Riemannian manifold (Al", g) to be isometrically 
embeddable in M^*^"). In particular, a fundamental, longstanding open 
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problem is to characterize intrinsic metrics on a two-dimensional Rie- 
mannian manifold Ai^ which can be realized as isometric immersions 
into M.^ (cf. [6l [131 [m [T7] the references cited therein). Important 
results have been achieved for the embedding of surfaces with positive 
Gauss curvature which can be formulated as an elliptic boundary value 
problem (cf. [6]). For the case of surfaces of negative Gauss curvature 
where the underlying partial differential equations are hyperbolic, the 
complimentary problem would be an initial or initial-boundary value 
problem. When the Gauss curvature changes sign, the problem then 
becomes an initial-boundary value problem of mixed elliptic-hyperbolic 
type. Hong in [7] first proved that complete negatively curved surfaces 
can be isometrically immersed in if the Gauss curvature decays at 
certain rate in the time-like direction. In fact, a crucial lemma in Hong 
[7] (also see Lemma 10.2.9 in [6j) shows that, for such a decay rate 
of the negative Gauss curvature, there exists a unique global smooth, 
small solution forward in time for prescribed smooth, small initial data. 
We are interested in solving the corresponding problem for a class of 
large non-smooth initial data. 

In Chen-Slemrod-Wang [3], we have introduced a general approach, 
which combines a fluid dynamic formulation of balance laws with a 
compensated compactness framework, to deal with the isometric im- 
mersion problem in (even when the Gauss curvature changes sign). 
In Chen-Slemrod-Wang |2j, we have developed a vanishing viscosity 
method to establish the existence of a weak entropy solution to the 
transonic flow in gas dynamics past an obstacle such as an airfoil, via 
the method of compensated compactness ([HI US])- We have found in 
[3] that the idea of [2j for gas dynamics is useful for solving the isomet- 
ric embedding problem in differential geometry. In particular, in [5], we 
have formulated the isometric immersion problem for two-dimensional 
Riemannian manifolds in via solvability of the Gauss-Codazzi sys- 
tem, and have introduced a fluid dynamic formulation of balance laws 
for the Gauss-Codazzi system. Then we have formed a compensated 
compactness framework and present one of our main observations that 
this framework is a natural formulation to ensure the weak continuity of 
the Gauss-Codazzi system for approximate solutions, which yields the 
isometric realization of two-dimensional surfaces in M^. As a first ap- 
plication of this approach, we have focused on the isometric immersion 
problem of two-dimensional Riemannian manifolds with strictly neg- 
ative Gauss curvature. Since the local existence of smooth solutions 
follows from the standard hyperbolic theory, we are concerned with the 
global existence of solutions of the initial value problem with large ini- 
tial data. The metrics (gij) we study have special structures and forms 
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usually associated with the catenoid of revolution and the helicoid. For 
these cases, while Hong's theorem [7] applies to obtain the existence of 
a solution for small smooth initial data, our result yields a large-data 
existence theorem for a C^'^ isometric immersion. To achieve this, we 
have introduced a vanishing viscosity method depending on the fea- 
tures of the initial value problem for isometric immersions and have 
presented a technique to make the apriori estimates including the L°° 
control and if "^-compactness for the viscous approximate solutions. 
This yields the weak convergence of the vanishing viscosity approxi- 
mate solutions and the weak continuity of the Gauss-Codazzi system 
for the approximate solutions, hence the existence of a C^'^-isometric 
immersion of the manifold into with prescribed initial conditions. 

From Chen- Slemrod- Wang [21 [3] , we have seen a remarkable connec- 
tion between the two distinct areas of gas dynamics and differential 
geometry. Here we present such a connection and show how the fluid 
dynamics can be used to formulate a geometry problem. Thus, we 
will present first the equations in Chen-Slemrod-Wang [2] for the tran- 
sonic flow problem in gas dynamics, and then the formulation using 
the fluid dynamic variables in conservation laws of gas dynamics for 
the isometric embedding problem in differential geometry. 



In two space dimensions with variables {x,y), the steady transonic 
flow of isentropic case is governed by the following steady Euler equa- 
tions on conservations of mass and momentum in gas dynamics: 



where p is the density, {u,v) is the velocity, and p = ^ {•y > 1) is the 
pressure. If we assume that the flow is irrotational, then system (12. ip 
can be reduced to the following two equations of irrotationality and 
conservation of mass: 



and, by scaling, the density p is determined by Bernoulli's law: 



2. Equations of Gas Dynamics 




(2.1) 






(2.3) 
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where q is the flow speed defined by = u'^ + v'^. The sound speed c 
is defined as 



c'=p\p) = l-^q'. (2.4) 
At the cavitation point p = 0, 



V 7 - 1 

At the stagnation point q = the density reaches its maximum p = 1. 
Bernoulh's law ( I2.3P is valid for < g < qcav At the sonic point q = c, 
fl2.4l) implies = Define the critical speed qcr as 



7 + 1 

We rewrite Bernoulli's law ( 12. 3p in the form 

- d = ^ - • (2-5) 

7 + 1 

Thus the fiow is subsonic when q < qcr, sonic when q = qcr, and 
supersonic when q > qcr- For the isothermal fiow (7 = 1), p = c^p 

where c > is the constant sound speed, the density p is given by 
Bernoulli's law: 

p = po exp ( ^-^) (2.6) 

for some constant po > 0, and qcr = c. 

3. Isometric Embedding in Differential Geometry 

In this section, we discuss the isometric embedding problem in dif- 
ferential geometry in and its formulation of fiuid dynamics. 

We first give the Gauss-Codazzi system of isometric embedding in 
W. Let gij,i,j — 1,2, be the given metric of a two-dimensional Rie- 
mannian manifold A4 parameterized on an open set O C M^. The first 
fundamental form / for on is 

I := gii{dxf + 2gi2dxdy + g22{dyf, 

and the isometric embedding problem is to seek a map r : — )■ 
such that dv ■ dv = /, that is, 

d^r ■ d^r = gu, d^r ■ dyV = g^, dyV ■ dyV = g22, 

so that {dxT, dyr} in are linearly independent. The corresponding 
second fundamental form is 

// := hii{dxY + 2hi2dxdy + h22idyY. 
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The fundamental theorem of surface theory (cf. \^ £j) indicates that 
there exists a surface in whose first and second fundamental forms 
are I and // if the coefficients (gij) and (hij) of the two given quadratic 
forms / and // with (gij) > satisfy the Gauss-Codazzi system. It is 
indicated in Mardare [lOj (Theorem 9; also see ^) that this theorem 
holds even when (hij) is only in for given (gij) in C^'^, for which the 
immersion surface is C^'^. This shows that, for the realization of a two- 
dimensional Riemannian manifold in with given metric (gij) > 0, it 
suffices to solve (hij) G L°° determined by the Gauss-Codazzi system 
to recover r a posteriori. The Gauss-Codazzi system (cf. [HE]) can be 
written as 

'd^M - dyL = T^^L - 2V'§M + T^^^N, 



d^N - dyM = -r^2^L + 2rS2^M - T^^^N, 



(3.1) 



with 
where 



/ N -^1212 J-, f T^im) o -r^im) . -r^{n)-r^{m) -p, (n) -p, (m) 

i^i^^y) = -T-^^ Rijki = 9im [dk^j - 9jT\^ + ^j'Kk - ^ik K/ 



LN -M^ = K, (3.2) 
r i\/r Ar ^22 I I , , . X 2 

" VM' " VH' " VM' " ^^^'^^^J''' = ^11^22 - gi2, 

k{x, y) is the Gauss curvature that is determined by the relation: 

\9\ 

Rijki is the curvature tensor and depends on {gij) and its first and 
second derivatives, and 

r!? = {d,ga + d,g,i - dig,,) 

is the Christoffel symbol and depends on the first derivatives of {gij)-, 
where the summation convention is used, {g^^) denotes the inverse of 
{gij), and {di, 82) = {d^, dy). Therefore, given a positive definite metric 
{gij) G C^'^, the Gauss-Codazzi system gives us three equations for the 
three unknowns {L, M, N) determining the second fundamental form 
II. Note that, although {gij) is positive definite, -R1212 may change sign 
and so does the Gauss curvature n. Thus, the Gauss-Codazzi system 
f l3.ip - fl3.2p generically is of mixed hyperbohc-elhptic type, as in tran- 
sonic fiow (cf. [2]). In Chen-Slemrod-Wang |3], we have introduced a 
general approach to deal with the isometric immersion problem involv- 
ing nonlinear partial differential equations of mixed hyperbolic-elliptic 
type by combining a fluid dynamic formulation of balance laws with a 
compensated compactness framework. As an example of direct applica- 
tions of this approach, we have shown how this approach can be applied 
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to establish an isometric immersion of a two-dimensional Riemannian 
manifold with negative Gauss curvature in R^. 

We now describe the fluid dynamic formulation of the Gauss-Codazzi 
system fl3.ip - fl3.2p in detail. Although, from the viewpoint of geom- 
etry, the constraint condition (13. 2p is a Monge-Ampere equation and 
the equations in (13. ip are integrability relations, we can put the prob- 
lem into a fluid dynamic formulation so that the isometric immersion 
problem may be solved via the approaches for transonic flows of fluid 
dynamics in Chen-Slemrod-Wang [2j. To do this, we set 

L = pv^ + p, M = —puv, N = pu^ + p, 

and set = + as usual. Then the equations in (13. ip become the 
familiar balance laws of momentum: 

' d,{puv) + dyi^pv" + p) = -(pt;2 + p)rg - IpuvV^^l - {pu^ + p)Tf^, 



d-Apu'+p) + dy{puv) = -(pt;2 +p)rg - 2puvT'^S " (pw' +p)r!l\ 

(3.3) 



and the Monge-Ampere constraint 03. 2p becomes 

ppq^ + = K. (3.4) 

We choose pressure p as for the Chaplygin-type gas: 

1 

p = . 

P 

Then, from (13. 4p . we have the "Bernoulli" relation: 

P = ^4=- (3.5) 

+ 1^ 

This yields 

p = -\/q^ + (3.6) 

and the formulas for and f ^: 

=p{p- M), v^ = p{p-L), = {N - p){L - p). 

The last relation for gives the relation for p in terms of (L, M, A^), 
and then the first two give the relations for {u, v) in terms of (L, M, A^). 
We rewrite (13. 3 p as 



dxipuv) + dy{pv^ +p) = Ri, 
dxipu'^ +p) + dy{puv) = R2, 



(3.7) 
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where -Ri and R2 denote the right-hand sides of fl3.3p . Then we can 
write down our "rotationahty-continuity equations" as 



In summary, the Gauss-Codazzi system fl3.1l) - fl3.2p . the momentum 
equations fl3.3p - fl3.6l) . and the rotationahty-continuity equations (13.50 
and (13. 8p are ah formahy equivalent. However, for weak solutions, 
we know from our experience with gas dynamics that this equivalence 
breaks down. In Chen-Dafermos-Slemrod-Wang p!], the decision has 
been made (as is standard in gas dynamics) to solve the rotationahty- 
continuity equations and view the momentum equations as "entropy" 
equalities which may become inequalities for weak solutions. In ge- 
ometry, this situation is just the reverse. It is the Gauss-Codazzi sys- 
tem that must be solved exactly, and hence the rotationahty-continuity 
equations will become "entropy" inequalities for weak solutions. 
We define the "sound" speed as: 



Hence, under this formulation, 

(i) when k > 0, the "flow" is subsonic, i.e., q < c, and system (13. 3p - 
([S3D is elliptic; 

(ii) when k < 0, the "flow" is supersonic, i.e., q > c, and system 
([S3D-(IS3D is hyperbolic; 

(iii) when k = 0, the "flow" is sonic, i.e., q = c, and system (13.30 - 
(13. 4p is degenerate. 

In general, system (I3.3p - (l3.4p is of mixed hyperbolic-elliptic type. 
Thus, the isometric immersion problem involves the existence of so- 
lutions to nonlinear partial differential equations of mixed hyperbolic- 
elliptic type. 

In Chen-Slemrod-Wang [5], we have considered one of the spatial 
variables x and y as time-like, have introduced a vanishing viscosity 
method via parabolic regularization to obtain the uniform L°° esti- 
mate by identifying invariant regions for the approximate solutions, 
and have shown that the i?;~^-compactness can be achieved for the 
viscous approximate solutions. Then, as in Chen-Slemrod-Wang [2|, 
the compensated compactness framework yields a weak solution to the 
initial value problem of system (I3.3p - (l3.4p when the initial data lies 




(3.8) 



c2 = p'{p) = i^, 



(3.9) 



then from our "Bernoulli" relation (13. 5p . we see 

2 2 I 
c = q + K. 



(3.10) 



8 



GUI-QIANG CHEN, MARSHALL SLEMROD, AND DEHUA WANG 



in the diamond-shaped invariant region. This estabhshes a C^'^(M^) 
immersion of the Riemannian manifold into M^. In particular, our ex- 
istence result asserts the existence of a C^'^-surface for the associated 
metric for a class of non-circular cross-sections prescribed at x = 
for catenoid. Our study in [3] also applies to the helicoid. See [3] 
for the details. Possible implication of our approach may be in exis- 
tence theorems for equilibrium configurations of a catenoidal shell as 
detailed in Vaziri-Mahedevan [TB]. However, the existence of isomet- 
ric embeddings/ immersions of a general surface with negative Gauss 
curvature is still open. When the Gauss curvature k changes sign, the 
problem becomes transonic and thus mixed hyperbolic-elliptic type. In 
this mixed-type problem, only special local solutions are known to ex- 
ist for special data ([316]), and the existence of global solutions is a 
significantly difficult open problem. 
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